arXiv:l 502.06749vl [math-ph] 24 Feb 2015 


One-dimensional two-component Bose gas 
and the algebraic Bethe ansatz 


N. A. SlavnovQ 


Steklov Mathematical Institute, Moscow, Russia 


Abstract 

We apply the nested algebraic Bethe ansatz to a model of one-dimensional two-compo¬ 
nent Bose gas with 5-function repulsive interaction. Using a lattice approximation of the 
T-operator we find Bethe vectors of the model in the continuous limit. We also obtain a 
series representation for the monodromy matrix of the model in terms of Bose fields. This 
representation allows us to study an asymptotic expansion of the monodromy matrix over 
the spectral parameter. 


1 Introduction 

In this paper we consider a model of one-dimensional two-component Bose gas with 5-function 
repulsive interaction (TCBG model). This model is a generalization of the Lieb-Liniger model 
0,0 (Quantum nonlinear Schrodinger equation), in which Bose fields have two internal degrees 
of freedom (colors). This model was solved by C. N. Yang 0 where the eigenvectors and the 
spectrum of the Hamiltonian were found. The general approach to the solution of the model 
with n internal degrees of freedom (multi-component Bose gas) was given in 0] (see also 0, 6;j]). 
The nested algebraic Bethe ansatz was applied to this model in [3,; 8|] . The main goal of this 
paper is to create a base for calculating from factors of local operators in this model in the 
framework of the nested algebraic Bethe ansatz. 

The algebraic Bethe ansatz is an efficient method for finding the spectra of quantum Hamil¬ 
tonians. However, in a viewpoint of calculating form factors of local operators application of 
this method meets some difficulties. The main problem is to embed the local operators of the 
model under consideration into_the algebra of the monodromy matrix entries. In some cases, 
this problem can be solved 00. However, to construct such a solution requires that the mon¬ 
odromy matrix of the model T(u) would be expressed in terms of the i?-nratrix. This is not the 
case of the TCBG model. On the other hand, in the framework of the traditional approach one 
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can easily obtain representations for form factors of local operators and correlation functions 
in terms of multiple integrals of the product of the wave functions. However, the evaluation of 
those multiple integrals is facing serious technical difficulties, and still they have been computed 
only for some relatively simple special cases 


Hi- 


Recently a method of calculating form factors of local operators in models possessing GL{ 3) 
symmetry was developed in [12]. This method is based on the nested algebraic Bethe ansatz 
and deals with partial zero modes of the monodromy matrix entries Tij(u ) [l3[ in a composite 
model 14]. Most of the tools of this approach can be directly used in TCBG model, however 
some of them should be slightly modified. In particular, one should adjust a definition of the 
zero modes. We solve these problems in the present paper. 

We consider a lattice approximation of the TCBG model. Using the T-operator obtained in 
pi IS] we construct a monodromy matrix and Bethe vectors. We show that these vectors have 
a correct continuous limit. We also obtain an explicit series representation for the monodromy 
matrix in terms of local Bose fields. Using this representation we are able to derive an asymptotic 
expansion of the monodromy matrix over the spectral parameter. In this way we find the zero 
modes. 

The paper is organized as follows. In section [2] we describe a general scheme of the algebraic 
Bethe ansatz. We define Bethe vectors of GL(3)-invariant models and give their representation 
in a multi-composite model. Section [3] is devoted to a brief description the TCBG model. In 
section 3] we give a lattice approximation of the TCBG model in the framework of the nested 
algebraic Bethe ansatz. In section [5] we consider continuous limit of the Bethe vectors of the 
lattice model. In section [6] we obtain a series representation of the TCBG monodromy matrix. 
Using this representation we find an antimorphism between Bose fields in section [7] and zero 
modes of the monodromy matrix entries in section [8] In conclusion we discuss some further 
applications of the results obtained. 


2 Algebraic Bethe ansatz 


In this section we describe an abstract scheme of the algebraic Bethe ansatz, which is valid for a 
wide class of quantum integrable models li|- 17]. The key objects of the algebraic Bethe ansatz 
are a monodromy matrix and R-matrix. The models considered below are described by the 
GL(3)-invariant R-matrix 18, 0 acting in the tensor product Ui <g) V 2 of two auxiliary spaces 
U fc -C 3 , As = 1,2: 


R(x,y) =1 +g(x,y)P, g{x,y) = 


x-y 


( 2 . 1 ) 


In the above definition, I is the identity matrix in Vj V 2 , P is the permutation matrix that 

exchanges V\ and V 2 , and c is a constant. 

The monodromy matrix T(w ) satisfies the algebra 


Rl2{wi,W2)Ti(wi)T 2 (w2) = T 2 {w2)Ti(w 1 )Ri2(wi,W2). ( 2 . 2 ) 


Equation (12.21) holds in the tensor product Ui (g> V 2 <S> H, where 7i is the Hilbert space of the 
Hamiltonian of the model under consideration. The matrices Tfc(rc) act non-trivially in Vk<S>T~L- 
We assume that the space % possesses a pseudovacuum vector |0). Similarly the dual space 
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7i* possesses a dual pseudovacuum vector (0|. These vectors are annihilated by the operators 
Tij(w), where i > j for |0) and i < j for (0|. At the same time both vectors are eigenvectors of 
the diagonal entries of the monodromy matrix 


Tn(w)\0) = Aj(u>)|0), = Ai(if;)(0|, * = 1,2,3, 


(2.3) 


where A i(w) are some scalar functions. In the framework of the general scheme of the alge¬ 
braic Bethe ansatz Xj(w) remain free functional parameters. Actually, it is always possible to 
normalize the monodromy matrix T(w ) —>• A^ 1 (rc)T(u;) so as to deal only with the ratios 


n(w) 


Ai jw) 
A 2{w) ’ 


r 3 (w) 


X 3 (w) 
X 2 (w )' 


(2.4) 


Below we assume that A 2 (tc) = 1. 

The trace in the auxiliary space V ~ C 3 of the monodromy matrix tr T(w) is called the 
transfer matrix. It is a generating functional of the Hamiltonian and all integrals of motion of 
the model. 


2.1 Notation 


We use the same notation and conventions as in the papers 0 , 0 - Besides the function g(x,y) 
we also introduce a function f(x,y) 


f(x,y) = l + g(x,y) = - —(2.5) 

x-y 

We denote sets of variables by bar: w, u, v etc. Individual elements of the sets are denoted 
by subscripts: Wj , Uk etc. Notation Hi , means u\ui etc. We also consider partitions of sets 
into disjoint subsets and denote them by symbol =>. Subsets are denoted by superscripts in 
parenthesis: u^\ For example, the notation u =>• u^} means that the set u is divided 

into two disjoint subsets and u^ 2 \ such that n u ^ = 0 and {u^ 1 \u < ' 2 ' > } = u. 

In order to avoid too cumbersome formulas we use a shorthand notation for products of 
operators or functions depending on one or two variables. Namely, if the operators T t] or the 
functions (12.4p depend on sets of variables, this means that one should take the product over 
the corresponding set. For example, 


T ij(u) = II TijiukY, r 3 (u (1) ) = r 3 (uj). (2.6) 

ui~£u ujeuM 

Similar convention is applied to the products of the functions f(x,y ): 

f(z,Wi)= f(z,wj ); f(u,v)=Y[ UfM ( 2 - 7 ) 

Wj €w Uj £u Vk Gh 

Wj^Wi 
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2.2 Bethe vectors 


The eigenvectors of the transfer matrix are called on-shell Bethe vectors (or simply on-shell vec¬ 
tors). In order to find them one should first construct generic Bethe vectors. In the framework 
of the algebraic Bethe ansatz generic Bethe vectors are polynomials in operators with i < j 
applied to the pseudovacuum vector. We denote them by M at b(u;v), stressing that they are 
parameterized by two sets of complex parameters u = {«i,..., u a } and v = {m,..., Vb} with 


a, b = 0,1,.... Different representations for Bethe vectors were found in 
one of the representations obtained in 


22-1251]. We give here 


B a ,b(u;v) = ^2 K "^!- ^ —^/(v (2) ,'c (1) )/(w (1) ,M (2) )Ti 3 (u (1) )T 23 (u (2) )Ti 2 (ii (2) )| 0 ). (2.8) 

Here the sums are taken over partitions of the sets u => u^} and v => (th 1 ),?^ 2 )} with 

0 < #u^ 1) = = n < min(a, b). We recall that the notation Ti 3 (u^) (and similar ones) 

means the product of the operators T\ 3 (u) with respect to the subset Finally, K n (i;( 1 ) l^ 1 )) 
is the the partition function of the six-vertex model with domain wall boundary conditions [26]. 
Its explicit representation was found in [2?| 


Kn(*| y) 


g(xj,xk)g(yk,yj) 


1 <k<j<n 


f(x,y ) / g 2 (xj,y k ) \ 

g(x,y) n \f( Xj ,y k )) 


(2.9) 


In particular, Ki(x|y) = g(x,y). 

A generic Bethe vector becomes on-shell, if the parameters u and v satisfy a system of Bethe 
equations: 


r\{ui) 

r 3 ( Vj ) 


f(Uj,Uj) 
f{Ui,Ui ) 


f(v,Ui ), 
f(vj,u), 


i — 1,..., a, 
j = 1,... ,6. 


( 2 . 10 ) 


Recall that Ui = u\ui and Vj = v \ Vj. 


2.3 Multi-composite model 

Study of the properties of local operators in the framework of the algebraic Bethe ansatz can 
be done by the use of a composite model 0- Suppose that we have a lattice quantum model 
of N sites. Then the monodromy matrix T(u ) is a product of local L-operators 

T[u) = Ln{u) ... L\(u). (2-11) 

Let us fix an arbitrary site m, 1 < m < N. Then (12.111) can be written as 

T(u) = T (2) (u)T (1) (u), (2.12) 

where 

T^{u) = L m {u)...L 1 {u), T^{u) = L N {u)...L m+1 {u). (2.13) 
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Representation (12.1211 defines a composite model. In the framework of the composite model 
the original matrix T(u ) is called the total monodromy matrix, while the matrices T^ 2 \u) and 
Tdl(u) are called partial monodromy matrices. The matrix elements of the partial monodromy 
matrices T^\u) and T^ 2 \u) act in the spaces and associated to the lattice intervals 
fl.ml and \m + l.JVl respectively. The entries of the total monodromy matrix act in the space 
of states « = 

In the framework of the algebraic Bethe ansatz it is assumed that and T-L ^ possess 
pseudovacuum vectors |0)( fc \ k = 1,2, such that |0) = |0) ( ' 1 ) <8> |0)^. These vectors have the 
properties analogous to (12.31) 


TjP(u)\0)W=0, i>j , 


T^(u)\0)^ = X\ k \u)\0)W, A: = 1,2. 


- \ (fc) 


(2.14) 


Similarly to (12.41) we introduce ratios 

\[ k \w) 


(fc) ( \ 
r 1 (w) = 


„(*) 


\ ( 2 k \w) 


H = 


x { 2 k \w) 


k = 1 , 2 . 


(2.15) 


Due to the normalization A 2 (u) = 1 we can always set X^ 2 \u) = 1. Below we also extend 
convention (12.61) to the products of functions (12.151) . 

One can construct for every partial monodromy matrix T^ k \u) the corresponding partial 
Bethe vectors B^('u: v). They are given by equation (12.81) . where one should replace all Tij{u) 
by T^ k \u) and |0) by |0)( fc ). The main problem considered in the framework of the composite 

model is to express total Bethe vectors B a ^(u]v) in terms of partial B^('u; v). This problem 
was solved in 114] for GL(2)-based models. More general case of GL(N) - invar i ant models was 
considered in |22.fi^|. Particular case of GL(3)-invariant models was studied in [j^ |. where the 
following representation was found: 


ia,b(u-v) = ^rf ) (/ ) ) 4 1 ) (i; (2) ) /(,J i(^ ( 1 ) ^ ( 1 ) ) B i 2 2 ) , 6 2 (« (2); ^ (2) 


(2.16) 


Here the sum is taken over all possible partitions u =>• u^} and v => }. The 

cardinalities of the subsets are shown by the subscripts of the partial Bethe vectors. 

Similarly we can define a multi-composite model, where the original interval is divided into 
M > 2 intervals 

T[u) = T^ M \u)...T {l \u). 


(2.17) 


For each of these intervals we can define partial Bethe vectors . Then the total Bethe 
vector can be expressed in terms of the partial ones as follows 


^a,b{u-v) = Y^ n 

1 <k<j<M 




f(yk>\ uW) 


M 


3 =1 

(2.18) 
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Here the functions r^\u) and r^\v) are vacuum eigenvalues of the operators T ^ 1 (u) and 
T.j.j'* (v) respectively. The sum in (12.181) is taken over all possible partitions 


Si) 


(j), 


u 


#u^=aj, a\ H-h ajf = a, 


# v ti) = bj, bi H-1 -b M = b. 


(2.19) 


It is important that the number M of the partial monodromy matrices is not related to the 
cardinalities of the Bethe parameters a and b. In particular, we can have M > a and M > b. 
In this case some of numbers a,j and bj are equal to zero, that is, the corresponding subsets are 
empty. 

Equation (|2.18l) can be easily proved by induction over M. Indeed, assuming that it is 
valid for M — 1 partial monodromy matrices we apply (|2.16l) to the partial Bethe vector 
]'l M tK M-1 )). This immediately gives (|2.18|> for M partial monodromy matrices. 

In particular cases a = 0 or b = 0 we reproduce known formulas for Bethe vectors in 
GL(2)-invariant multi-composite model 0, [in| For instance, 


M 


1,0 = 


»(“)=e n 

l<k<j<M 


r[ j) (u^)f(u® , u {k) } Yl IB g) ( u& ). 

j =i 


( 2 . 20 ) 


The multi-composite model is a convenient way to express the Bethe vectors in terms of local 
operators. In the next section, we discuss the method in more detail. 


2.4 Bethe vectors in the SU( 2) XXX chain 

As a first application of the multi-composite model we construct Bethe vectors of the SU( 2) 
inhomogeneous XXX chain. This result will be used in section [3] for description of Bethe 
vectors of TCBG model. 

Consider an inhomogeneous XXX chain consisting of M sites. This model has a 2 x 2 
monodromy matrix T- xxx ) (^). therefore Bethe vectors are parameterized by only one set of the 
Bethe parameters, say u. Respectively, considering the multi-composite model one should use 

firm 

The monodromy matrix is defined as a product of local L-operators 

T^ xxx \u) = L^ xx) (u - £ m ) • • • L {xxx \u - 6), (2.21) 


where are inhomogeneities and 



I A*+§(1 + 0 ca n A 

u\ ca+ «+§(! -<)/ 


( 2 . 22 ) 


Here and are spin-1/2 operators acting in the n-th site of the chain. They are given 
by the standard Pauli matrices acting in the n-th copy of the tensor product (C 2 ) " . The 
pseudovacuum vector is the state with all spins up 


10 ) — (o) Af ® ® (o) l • 


(2.23) 
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Bethe vectors with a spins down and M — a spins up have the form 



E n-* 

M>j a >—>j i>l m=1 


a- 


| 0 ). 


(2.24) 


where (it; £) are coefficients depending on the Bethe parameters u and inhomogeneities 

£. Let us find these coefficients explicitly. 

Consider a multi-composite model with M partial monodromy matrices It means that 
every coincides with the L-operator Lj[u — £j). Then every partial Bethe vectors B (u^) 
in (12.201) corresponds to the j-th site of the chain, therefore due to (12.221) we obtain 

= g(u(o-y (1).. (2.25) 

Obviously vanishes if aj > 1, because (a)”) 2 = 0. Thus, we conclude that aj < 1 and the 
subsets are either empty or they consist of one element. Let subsets (k = 1,... , a) 
corresponding to the lattice sites ji,...,j a contain one element u k , while other subsets are 
empty. Then the sum over partitions of the set u turns into the sum over permutations in u 
and the sum over the lattice sites ji, ■ ■ ■, j a with the restriction j a > ■ ■ ■ > ji¬ 
lt is easy to see that 


u ~ £j + I(1 + a j) 
u-tj 




u - £j + §(! - 

u-Zj 


a)i = a\ 


(2.26) 


and thus, 

r[ j) (u) = (2.27) 

Then equation (12.201) takes the form 


®<r*)(u) = Sym [] f( Uj , 

u k) n 

M 

( n 

f(Uk,£m))9(u k ,t jk )(T jk 

ZL . ., , . 

1 <k<j<a 

M>j a >- >ji>l k= 1 

m=j k + 1 



(2.28) 

where the symbol Sym means symmetrization (i.e. the sum over permutations) over the set 
indicated by the subscript. The symmetrization in (12.281) acts on all the expression depending 
on u. Comparing (12.281) with (12.241) we see that 


’, M ] a («; 0 = s y_ m II /(“i> u k) n 


1 <k<j<a 


k =1 


M 


(n 

f ( u ki ^m)^j S( u kt Cjk ) 

m=j k +l 



(2.29) 


In the homogeneous limit = c/2 this expression coincides with the amplitude of the Bethe 
vector in the coordinate Bethe ansatz representation (see 0 )- 
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3 Two-component Bose gas 

We consider the TCBG model on a finite interval [0, L\ with periodic boundary conditions. In 
the second quantized form the Hamiltonian has the form 

H = £ (d x ^ a d x ^ a + dx, (3.1) 

where x > 0 is a coupling constant, ck,/ 3 = 1,2 and the summation over repeated subscripts is 
assumed. Bose fields T^x) and 'Fa(x) satisfy canonical commutation relations 

^(y)] = 5a/3$(x - y)- (3.2) 

The coupling constant x is related to the constant c in (12.111 by x = ic. 

The basis in the Fock space of the model is constructed by acting with operators tJ^x) onto 
the Fock vacuum |0) defined as 

*«(x) |0)=0, <0|Tt (x) = 0, (0|0) = 1. (3.3) 

Observe that in the case of the TCBG model the pseudovacuum vector (12.311 coincides with the 
Fock vacuum 10), therefore we use the same notation for them. 

The spectral problem for the TCBG model was solved in Q (see also (4, 6]). The Hamiltonian 
eigenvectors can be found in two steps. Using the terminology of the algebraic Bethe ansatz 
one can say that at the first stage one should construct a generic Bethe vector B v). In the 
TCBG model Bethe vectors exist for a <b. They have the following forirj^]: 


,b(u;v)= ^ J dzi...dz b Xk 1 ,...,k a (zi,---,z b \u,v ) 


b>ka>--->kl>l 


V 


n n ®k^)io). (3.4) 


m= 1 


1=1 


Here the integration domain is T> = L > Zb > ■ ■ ■ > z\ > 0. In this domain the wave function 
Xk 1 ,...,k a (z 1 ,...,z b \u,v) has the form 


Xk u ...,k a (zi,...,Zb\u,v) = Symfi^jjii^ + c) f(vj,v k )Y[ 1 

V b>j>k> 1 k =1 


J'Z'kVk 


(3.5) 


where the coefficients (v/,v + c) are given by (12.2911 . 

Generic Bethe vector (13.411 becomes an eigenvector of the Hamiltonian (13.111 if the parameters 
u and v satisfy the system of Bethe equations (12.1011 . In the TCBG model it has the following 

2 Here and below we do not take care about eigenvectors normalization in all formulas for them. 
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form 0]: 


e iLvj _ 


n 

k=\ 

Ml 


Vj — V k + ix -|-r U£ ~ Vj + ix 


-n 


k=1 V 3 -V k -lX “ - Vj 


i=n 


a Ui — U£ — ix -r^r v k — Ui — ix 


n 


f=i Ui - Ul + IX “ v k - Ui 

Ml 


j = l,...,b, 


i = 1,..., a. 


(3.6) 


Comparing this system with (I2.10|) we conclude that in the TCBG model r\(u) = 1 and r%(u) = 
ALu 


4 Lattice two-component Bose gas 

Quantum systems describing by the GL(3)-invariant 77-matrix (12.11) were considered in 
There a prototype of a lattice L-operator of the TCBG model was found. It has the following 
form: 

L,( a \u) = ul +p, (4.1) 


where 


ajai 


aja 2 ia\ \Jm + p\ 

V = | a^ai a\a 2 ia\y/m + p I • (4-2) 

\iy/m + p a\ iy/m + p a 2 — m — p ) 

Here m is an arbitrary complex number and p = a\ai + 0 ^ 02 - The operators a k and a' k (k = 1,2) 
act in a Fock space with the Fock vacuum |0): afc|0) = 0. They have standard commutation 
relations of the Heisenberg algebra [aj,ajj = 5i k . 

The L-operator (14.11) satisfies the algebra (12.21) with R- matrix (12.11) at c = —1. Basing on 
the L-operator (14.11) one can construct a quantum system of discrete bosons. In order to obtain 
a continuous quantum system one should make several transforms of SI]). First, we introduce 
operators 

ipk = A ~ 1/2 a k , ipl = A -1/2 (4, k = 1,2, (4.3) 

so that 


h All = -£■ 


A 


(4.4) 


In these formulas A is a lattice interval. Setting m = we introduce a new L-operator as 

+ ^ (4.5) 

where J = diag(l, 1,-1). Obviously L(u ) satisfies the LTT-relation (12.21) with L-matrix (12.11) 
at c = —ix. 

The last transformation is to make N copies L n (n = 1,...,N) of L-operator (14.51) by 
changing ip k —> ip k (n) and ip' k —> i^ k (n) with 

SjkSnrn 




A 


(4.6) 
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The operators ipk(n) and V4( n ) are lattice approximations of the Bose fields 'Ffc(x) and T|,(x). 
Indeed, let us divide the interval [0, L] into N sites of the length A. Setting x n = nA and 


1 f Xn 

= — / ^ k (x)dx, 
^ J a: n _i 



1 

A 



(4.7) 


we reproduce commutation relations (14.61) . On the other hand, in the limit A —)■ 0 the operators 
(14.71) obviously turn into the Bose field^f T^x) an cl T^x). 

Now we can define a monodromy matrix in a standard way 


T(u) = L n {u) ... Li(u), 


(4.8) 


where 


L n {u) — 




and 


1 - 


iAQn'i/ji (n) 


^ + i£ y^V4( n )V ; 2( ri ') 

iAQn^in) 


—iAijj\(n)Q r 

—iA'll^ipjQr. 


f + H^ + 


^~Pnj 


(4.9) 


N= 1 - 


iu A 


p n = i\{n)^ 1 {n) + ^l{n)^ 2 {n), Q n =[^ + 


>c 


: A 2 


" Pn 


1/2 


(4.10) 


The normalization factor J\f in (14.91) is used in order to satisfy the condition A 2 (u) = 1. 

Remark. We write the number of the lattice site n as the argument of the operators and 
i/jj. Traditionally this number is written as subscript of i/'j and ipj, but in the case of the TCBG 
model it is not convenient. 

L-operator (14.91) is a natural generalization of a 2 x 2 L-operator found in [ 32 I ] for the lattice 
model of one-component bosons: 


L n {u) 


J_ /1 - ^ (n)i/j(n) -iA^{n)Q n 

AT iAQ n i/j(n) 1 + ^ (n)ip(n) 


(4.11) 


It is easy to see that L-operator (14.111) is the right-lower 2x2 minor of the matrix (14.91) with 
the identification 1^1 (n) = 0, 1^2 (ji) = ip(n). It was shown by different methods in [33M35f| that 
in the continuous limit A —>• 0 the L-operator (14.111) describes the model of one-dimensional 
bosons with 5-function interaction. We have to solve an analogous problem: to check that in 
the continuous limit the model with the monodromy matrix (14.81) and L-operator (14.91) does 
describe the TCBG model. For this purpose we will find Bethe vectors of the lattice model 
(14.81) and will show that they coincide with the states (13.51) in the continuous limit. 

Let us point out sever properties of the L-operator (14.91) . It is easy to see that 


(L n (u)) n |0) = (L n (w)) 22 |0) = |0>, (L„(u)) 33 |0) = r 0 (u) |0); 

i L n(u)) 12 |0) = 0, (0| (L n (tt)) 21 = 0, 

3 Here and below limits of operator-valued expressions should be understood in the weak sense. 
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where 


(l + mA\ 

ro(n)=[ rn ^). (4.13) 

From these properties we easily find 

ri(u ) = 1, 

Ti 2 (?x)|0) = 0, 

Note that in fact the condition r\(u) = 1 implies the actions of T\ 2 (u) and T 2 \ (u) in the second 
line of (14.141) . Indeed, we have from the RTT- relation (12.21) 

[T 21 (y),T 12 (u)\ = g(v,u)(T u (u)T 22 (v) - T u (v)T 22 (u)). (4.15) 

Applying this equation, for example, to the vector |0) and using r\(u ) = 1 we obtain 


rs(u) = (u); 

(0|T 2 i(u) = 0. 


(4.14) 


[T 21 (v),T 12 (u)}\0) = T 21 (v)T 12 (u) |0) 

= g(v,u)(T u (u)T 22 (v) -T u (v)T 22 (u))\0) = (n(u) - ri(u))|0) = 0. (4.16) 

Similarly, acting with (14.151) on (0| we obtain (0|T2i(w) = 0. 

The property Ti 2 (ii)| 0) = 0 leads to a simplification of the explicit formula for the Bethe 
vector (12.81) . Obviously, in this case we should consider only such partitions of the set u that 
= 0, and = u. Then (12.81) turns into 

M atb (u-,v) = Ka }(- f(^ 2 \v^)T 13 (v W )T 23 (v^)\0). (4.17) 

Here the sum is taken over partitions of the only one set v =4- u( 2 )} with a restriction 

= a. The last restriction evidently can be satisfied if and only if a < b. Hence, if a > b, 
then B ai fc(zZ; v) = 0. In particular, 

B o ,i(0;u) = T 23 (v) |0), Bi,i(«;u) = T 13 (v) |0). (4.18) 

To conclude this section we give two formulas concerning the continuous limit A —>• 0. The 
hrst formula gives the limit of powers of the function ?’o (u) 

/1 _l_ iuA \ Xn / A 

lim r!}(u) = lim . 2 = e iuXn . (4.19) 

A—5-0 0 V ’ A^0 ^ 1 - J K J 

The second formula describes a typical procedure of taking continuous limit of sums over the 
lattice sites. Let 4>(x) be an integrable function on the interval [0, L]. Then 


N N rXj 

3 =1 3 =1 Jx i- 1 


T((x) 


dx 



A 0, (4.20) 


and we recall that all limits of operator-valued expressions are understood in a weak sense. 
Thus we can formulate a general rule: a sum over the lattice sites multiplied by A turns into an 
integral in the continuous limit. It is easy to see that if we have an m-fold sum over the lattice 
sites multiplied by A m , then it turns into an m-fold integral in the continuous limit. 
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5 Bethe vectors in terms of local operators 


Consider a multi-composite model with the total monodromy matrix (|4.8|> . Let the number M 
of the partial monodromy matrices coincides with the number N of the lattice sites. Then every 
partial monodromy matrix T^ n \u) is the L-operator L n (u) m- Respectively we have 

r[ k \u) = 1, r^ k \v) = r 0 (v). (5.1) 

The formula for the total Bethe vector (12.181) takes the form 


N 


1)b (u;v) = 5^11 
3 = 1 


^ i 


(® 0) ) 




1 <k<j<N 




O ^0')). 


a j ,b 


(5.2) 


3 = 1 


This is the main formula that we shall use. But before applying this formula to the TCBG 
model it is useful to look how it works for a more simple example of the one-component Bose 
gas. 


5.1 One-component Bose gas 


The L-operator of the one-component Bose gas is given by (14.111) . however for the construction 
of Bethe vectors we need to know this L -operator only up to terms of order A: 


( 1 - 


_ iuA 

( ‘A /—‘ft \ 


-j i mA 
1 2 


+ 0(A 2 


(5.3) 


Recall that there we have set ip 2 {n) = ip(n ), = 0 and similarly for 'ifr,. (n). In the continuous 

limit these operators respectively turns into Bose fields 'F(x) and Tl(x). 

Bethe vectors of the one-component Bose gas correspond to the particular case of v) 

at a = 0 and u = 0. Then the formula (15.2p takes the form 


N N 

B o ,b(0;ti) = B b (v) = ^ JJb^u^). (5.4) 

j=1 1 <k<j<N j =1 

A partial Bethe vector in the site j is 

B^u^) = (-iA^i/j\j)) bj \0), (5.5) 

where corrections of the order 0(AL +1 ) are neglected. 

Remark. Recall that in the multi-composite model the total pseudovacuum vector |0) is 
equal to the tensor product of the partial pseudovacuum vectors |0)^ (j = 1,..., N). However, 
in the case of the one-component Bose gas we can assume that all operators ip(j) and act 

in the same Fock space. Obviously, due to commutativity of 'ip(j) and at j ^ k such the 

formulation is equivalent to the original one. In the case of the TCBG model we will use the 
same treatment of the multi-composite model. 

Consider an example 6 = 2. Then we have two possibilities. 


12 







• There exists one bj such that bj = 2, while all other bp = 0. Then the subset coincides 
with the original set {ui,u 2 }, while all other subsets jjW are empty. 

• There exist two bj and bk such that bj = bk = 1, while all other bi = 0. Then the subsets 

and v^ consist of one element (say, = V 2 and v^ = v\ or vice versa). All other 
subsets are empty. 

Consider the first case. We denote the corresponding contribution to the Bethe vector by 
1 B>2,0- Then 

N 

®2,0 = -^A 2 ^(r 0 (ui)ro(u 2 ))' ;_1 (V ,t (j)) 2 | 0 ) ) (5-6) 

3 =1 

and due to (I4.19|) we obtain 

N 

B 2i 0 = -xA 2 J2 e iXj(vi+V2) ^Hj)) 2 \0)- (5.7) 

3 = 1 

This sum goes to zero, because it has the coefficient A 2 . Indeed, due to (14.2011 we have 
N ,L 

A 2 'S~'e iXj ( Vl+V2 ' ) (ip\j)) 2 \0)—> A / e ix ^ +V2) (^(x)) 2 dx\0) ^ 0, A -> 0. (5.8) 

3=1 J ° 

It remains to consider the second case. We denote the corresponding contribution to the 
Bethe vector by Bu g. Then 

®i,i,0 = -^A 2 Sym rl\v2) r t\n)f( v 2,v 1 )^(j)^(m, (5-9) 

V l<k<j<N 

or due to (14.191) 

B ljl>0 = -xA 2 Sym £ e™ + ^ 2 /(u2 , Vl )^( k )\0). (5.10) 

v 1 <k<j<N 

This time we have again the coefficient A 2 , but the sum is double. Therefore the limit is 
finite 

rL r x 2 

lim B x i 0 = B 2 (u) = — xSym/(ii 2 , ui) / dx 2 / dxie lXlVl+1X2V2 T^(x 2 )T'l'(a:i)|0). (5.11) 

’ ’ v Jo Jo 

It is clear from (15.41) and (15.51) that for general b the Bethe vector B&(u) is proportional to A b . 
In the continuous limit this coefficient should be compensated. The only possible way to obtain 
such the compensation is to have a 6-fold sum over the lattice sites. Then A b times 6-fold sum 
gives a 6-fold integral. Hence, we should consider only such partitions of the set v = {ui,..., u&} 
that reduce to 6-fold sums over the lattice sites. Obviously, these are such partitions in which 
there are exactly 6 nonempty subsets. In this case every such subset consists of only one variable. 
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Thus, actually we deal with the case already considered in section l2~4l Therefore the sum over 
partitions reduces to the sum over the lattice sites and the sum over permutations, i.e. to the 
symmetrization over v. 

Thus, we obtain for general b 


N b 

= (-n/xA) b Sym f(vj,v k ) ^ JJ (r^ fc_1 (u fe )V’ t (jfc)) |°), (5.12) 

V b>j>k> 1 jb>—>ji k =1 

or partly taking continuous limit 

N b 

B&(u) = (-iv/^A) 6 Sym f{vj,v k ) ^ |0). (5.13) 

V b>j>k> 1 jb>" >ji k =1 

This 6 -fold sum over the lattice sites goes to a 6 -fold integral, and we finally arrive at 


lim BUu) 

A—>-0 



b 

d Xl ...dx b \[(e ix ^(x k j) 10 ), 

fc=l 


(5.14) 


where D = L>xj ) >--->x i > 0. Representation (15.141) coincides with well know result for 
the Bethe vectors in the coordinate Bethe ansatz Hal I 31 I ]. Thus, we have constructed Bethe 
vectors in terms of the local Bose field Tl(x) starting from the lattice L-operator (15.31) . 


5.2 Two-component Bose gas 

The infinitesimal lattice L-operator of the TCBG model has the form 0] 


L n {u) 


( 1 -^ 0 -iA^3hJ>\(n)\ 

0 1-^ -*Av^|(n) +0(A 2 ). 

yiAy/Te^iin) iAs/xip 2 {n) 1 + ^ J 


(5.15) 


We again consider a multi-composite model with the number of the partial monodromy matrices 

T (n) 

(u) equal to the number of the lattice sites. Then every T^ 1 ’ (u) of such the model coincides 
with the L-operator (15.151) . First of all let us find how Bethe vector depends on A. In the TCBG 
model Bethe vectors are given by (14.171) . It is easy to see that the total number of creation 
operators T 13 and T 23 in (14.171) is 6 . In the case of partial Bethe vectors b we 

have 

Ti 3 (w) = -iAy/iaf)\(j), T 23 (w) = -iAy/x^\{j). (5.16) 

Therefore 

^a-jbj ~ A bj , and thus, B aj 6 (u; u) ~ A 6 . (5-17) 

The Bethe vectors of the multi-composite TCBG model are given by (15.21) . Using the same 
arguments as in the case of the one-component Bose gas we conclude that we should consider 
only such partitions of the set v, where we have exactly 6 nonempty subsets consisting of one 
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element. Then the sum over such partitions of the set v turns into the sum over permutations 
of v and a 6-fold sum over the lattice sites. 

Consider now what happens with the partitions of the set u. In every partial Bethe vector 
bj > a,j. As we have shown above, all bj are equal either to zero or to one. If bj = 0, then 
a.j = 0. However if bj = 1, then either a,j = 1 or aj = 0. In the first case we obtain a partial 
Bethe vector of the form B x (, in the second case a partial Bethe vector of the form F?q {. But 
since all nonempty subsets consist of exactly one element, the sum over partitions of the set u 
also turns into the sum over permutations in u and a sum over the lattice sites where aj = 1. 

Thus, the sum in (15.21) is organised as follows. First, we should choose a set J consisting of 6 
numbers J = {ji,... ,jb}- These are the numbers of the lattice sites, where bj k = 1. In all other 
sites bj = 0. We assume that the subset v^ k ~) consists of one element ■ Taking symmetrization 
over v and the sum over all possible j k with the restriction jb > ■ ■ ■ > j l we thus reproduce 
the sum over partitions of the set v. More precisely, we reproduce only such partitions that 
eventually contribute into the continuous limit. 

Up to this point everything is exactly as in the case of one-component bosons. Now we 
should take into account the partitions of the set u. For this we should choose among the set 
J = {j\ ,..., jb} a subset of numbers K consisting of a elements: K = ,..., jk a }-> K C J. 

These are the numbers of the lattice sites, where a,-. = 1. In all other sites a ; = 0. We assume 

that the subset consists of one element u m . Taking symmetrization over u and the sum 

over all possible jk rn with the restriction j^ a > ■ ■ ■ > we reproduce the sum over partitions 
of the set u. 

Summarizing all above we recast (15.21) as follows: 


N 

B ai b(u; v) = Sym n f{Vj,V k ) JI f(uj,u k ) E E 

’ b>j>k> 1 a>j>k> 1 jb>—>31 ik a >—>3k-y 

jkm s J 

x n n f~ 1 ( v ^ u rn)]Jr J 0 k ~ 1 (v k )flM^\um]v krn ) JJ B^j ) (0;u £ ). (5.18) 

m= 1 e=k m +l k= 1 m= 1 je€J\K 

Due to (14.181) and (15.151) we find 

®o!i(0;u) = -iAy/m(>l(j)\0), B^|(u;u) = -iA^x ip\(j)\0), (5.19) 

J \ v i u ) 

and using ()4. 191) we obtain 

N 

K a ,b(u;v) = (-iA^) b Sym JJ f(vj,v k ) JJ f(uj,u k ) ^ ^ 

V,U b>j>k> 1 a>j>k> 1 jb>—>jl 3k a > —>jk 1 

jk m 

xfi ii n < 5 - 2 °) 

m=U= krn +l k =1 m=l j e eJ\K 
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Using obvious properties of the functions g(x,y) and f(x,y) 


f(x,y + c) = 


f(y,x )’ 


g(x,y + c) = - 


g{y,x) 
f(y,i 0 ’ 


(5.21) 


we see that 


9( v km j u n 


Sym l\ f( Uj ,u k ) JI ] . 

U a>j>k>l m.=lV J ^ Vk ™' Um > t=k m +l 

where the coefficients ka are given by (|2.29p . 

Thus (I5.20P takes the form 


n,/(^} =( - ir < 6, * (Sii+c) ' (5 ' 22) 


N 


*,b(u]v) = (-l) a (-iAy/x) b S ym f(vj,v k ) ^ ^ JJ e* 


IXkVk 


b>j>k>l 


jb>->ji 3ka>->3k 1 k =1 

jkm^J 


xn k a u!.,k a ( u 'i v + c ) n 4o»i°)> ( 5 - 23 ) 

m=l 


and it becomes clear that in the continuous limit we arrive at (ESP up to a normalization factor. 


6 Representation of the monodromy matrix in terms of Bose 
fields 

In this section we derive explicit representations of the monodromy matrix elements Tij{u ) in 
terms of the Bose fields. These representations have the form of a formal power series in the 
coupling constant x. It is worth mentioning, however, that in a weak sense, these series are cut 
on an arbitrary Bethe vector. 

Let us present the infinitesimal L-operator (I5.15P as a block-matrix of the size 2x2: 

Ln(u) = ( “ h 2)+0{ A 2 ). (6.1) 

Here cL = 1 + pp, and a is a 2 x 2 matrix a = (1 — pp) • 1, where 1 is the identity matrix of 
the size 2 x 2. A two-component vector-column b n and two-component vector-row Cn are 

b n = -iA^Jx () , Cn = iAi/x th(n)) . (6.2) 

\W n )/ 

It is convenient to separate the diagonal and anti-diagonal parts of the L-operator (16.111 as 
follow^: 

L n {u) = A (it) + W n , (6.3) 

4 Here and below we omit the terms 0(A 2 ) as they do not contribute to the continuous limit. 
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where 


A(u) = 


a 0 
0 d) ’ 


W n = 


0 b n 

"n 0 


(6.4) 


Now representation (16.31) should be substituted in (14.8|) and then developed into the series 
in W n . Since the anti-diagonal part W n is proportional to ^/x, the monodromy matrix T(u) 
becomes a polynomial in ^/x, that turns into an infinite power series in the continuous limit: 


T(u) = Y,x n/2 T n (u), 


(6.5) 


n=0 


where 


x n/2 T n (u) =(l- 


9 7 / A \ ~ N 

^ A N -k-W kri A kn ~ kn - 1 - 1 ■ ■ ■ ( 6 . 6 ) 


N>k„>->ki>l 


It is clear from (16.61) that the diagonal blocks of the monodromy matrix are series in integer 
powers of x, while anti-diagonal blocks are series in half-integer powers of x. 

Let 


— A _ ( 0 b ki 


W ki = A - Ki W ki A Ki ~ L = 


c ki 0 


where 


bki — 


h. 


iuA 


-I lULA l 1 IULA 

1 2 \ 1 2 . 


°ki = 


c ki 


1 - 


iuA 


1 IULA 1 1 I IULA 

1 2 2 , 


Then equation (16.61) takes the form 

*”/ 2 r„(u) = ( A a 


N>k n >--->ki>l 


(6.7) 


( 6 . 8 ) 


(6.9) 


Partly taking the continuous limit via (14.191) we obtain 

-N 


lim (l-iAV 

A^o \ 2 / 


A^ = 


0 e 


iuL i ’ 


and 


( 6 . 10 ) 

b k . = b ki e luXk i, c ki = c ki e~ luXk i. (6.11) 

It is convenient to study the operators T n (u ) separately for n even and n odd. Let n = 27. 
The product of two matrices W ki and W ki _ 1 gives a block-diagonal matrix 


W ki W ki _ i = 


_ / bk i c ki _ 1 


0 

Cki bki-1 


( 6 . 12 ) 


x e T 2l (u) = 


Ae(u) 0 

0 Dz(u) J 


(6.13) 
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Hence, we obtain 











where 


(6.14) 


and 


A((u) ] bk2e.Ck2t-1Pk2e.-2Ck2e.-3 ' ' 'bk 2 Cki j 

N'>k2£> • • ‘>k\ >1 


A?(«) = e*“ L 5^ c k 2eb k 2e-iC k 2e-2bk 2l - 3 ■'' Cfc 2 &ki- (6.15) 

N>k2£ >•••>&! >1 

Observe that all operators in these products commute, because they are from different lattice 
sites. Therefore 


bk 2e Ck 2 e_Pk 2 e_ 2 Ck 2 e_3 ’ ’ ’ bkPki — ■ bk 2e Ck 2e _Pk 2 e_ 2 c fc 2 ^_ 3 ''' bfaCfa ., 
Ck 2 pk 2 e_Pk 2 e_ 2 bk 2 e-3 ''' CkPk 1 — • Ck 2 eb k 2i-pk 2 (- 2 bk 2 e-3 ‘'' CkPk 1 ■> 
where the symbol : : means normal ordering. Obviously 

Ck 2 pk 2i -i = xA 2 e m(xfc 2i-i“ a: *2i) : (v 4 (& 2 i-i)V’ife) + Ip2p2i-i)^2{k2i)) ■ ■ 
Thus, we find 


(6.16) 


(6.17) 


APu) = ^A 2£ Ye 

iV>^>->/ci>l i=l 


x : mw(w*<*«>+*<«*<*-» ($£&$ $£&£>) : • (6 ' 18) 


2=1 


and 


D t (u) = e iuL y* A n Y n e ' m(Xfc2l ‘ Xfc2l -i) 

N>k 2 e> -->ki>l i =1 


X : 11(^1 fe-OVhfe) + V4( fc 2i-l)V , 2(fc2i)) : ■ (6-19) 


i=l 


It remains to replace the sums over ki by integrals via (14.201) . It is convenient to set Xk 2i = 2* 
and = y t . Then 

rL t 


Apu) = | e m(2i ^ etej c%} 0*(z, y) 


/0 )=i 

£-1 


x 


=n ( .!w.^o + .!„)..(*.«)) =. <«■»> 


and 


rL < 


£>*(«) =e iuL x^ / n{ e - m(2 ' l ‘ y,:) ^^} eK ^ : ^ : n(*!(^iW + *t( W )* aW ) : 

■1° i=i 1 i =1 


( 6 . 21 ) 
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where 


e -1 


@e(z, y ) = 0{zt - Ve) %*+1 - 2i)0(zi - y*). 


( 6 . 22 ) 


i=l 


The anti-diagonal blocks of the monodromy matrix can be found exactly in the same manner. 
Setting n = 2£ + 1 in (16.91) we find 


= ( c ° u) B i u) ) 


(6.23) 


where 


„L £ 

Ce(u) = ie mL x e+1/2 / dzi %} e ~ luye+1 dy e+1 d(y e+1 - zg)Ge(z, y) 

»=i L 


x: IK«K Zi)^i(y i+l ) + 'S>\(zi)'S> 2 {yi+i)) ■ (^i(yi); ^ 2 ( 2 / 1 )) : , (6.24) 


1=1 


and 


pL * 

Be{u ) = -ix £+1/2 / JJ | e -™(^-2/d dzt dyij e luye+1 dy e+ i9(y £+1 - ze)Gg(z, y) 

x : + 'S>l(y i )'S' 2 {zi)) • : • (6.25) 

Thus, we have obtained the explicit series representation for the monodromy matrix entries 
Tij{u ) in terms of the local Bose fields. This series is formal, and we do not study the problem 
of its convergence. It is easy to see, however, that if we introduce a vector 

l‘L a b 

|$a,&) = / dxi,..., dx a dyi,..., dy b $ a , b (xi,... ,x a \yi,... ,y b ) n»lwn*!wm. 

i =1 j =1 

(6.26) 

where & a ,b(x 1 , • • • ,x a ‘,yi ,... ,y&) is a continuous function within the integration domain, then 
the action of any Tij(u) on |$ a t) turns into a finite sum. 

7 Mapping of fields 

Due to the invariance of the 12-matrix under transposition with respect to both spaces, the 
mapping 

<j){T jk (u)) = T kj (u) (7.1) 

defines an antimorphism of the algebra (|2.2I) (see 0). This mapping is a very convenient tool 
in studying form factors, because it allows one to relate form factors of different operators. In 
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the case of the TCBG model antimorphism (17.11) agrees with the following mapping of the Bose 
fields: 

^(*i(s)) = ~*\( L - x), ^\(x)) = ~^i{L - x). (7.2) 

Indeed, consider, for example, how the mapping (17.21) acts on the matrix elements Tj k (u) for 
j,k = 1,2. Due to equations (16.131) . (16.201) we have 

OO 

Tjkiu ) = j/(T 2i ) jk (u), j, k = 1,2, (7.3) 

£=o 


where 

,L t t-1 2 

(T 2 e) jk {u) = / n{ e m ^- yn) dz n dy n }Q £ (z,y) :^ j (z £ )'& k (y l )Y[(^2'l>l(z n )'S> s (y n+ ij S ) : . 

4° n =1 k n=1s=l 

(7.4) 

Recall that due to the factor &i(z, y) the integral in (17.41) is taken over domain z^> yi > zp_\ > 
■ ■ ■ > zi > r/i. Therefore all the operators in (ED commute with each other, and actually the 
normal ordering is not necessary. Acting on (17.41) with <fi as in (17.21) we obtain 

pL (■ 

<K(T2t) .-»)= / \\{e iu ^-y^d Zn dy n \Q i {-z,y):^ k {L-y 1 )^ ] (L-z,) 

40 n=1 L J 

t -1 2 

*n (e Tt(L - y„ + i)T s (L - 2n )) : . (7.5) 

71—1 S = 1 

Now it is enough to change the integration variables z n —>• L — and —>• L — Z£+i-n- 

Then we have 


®e(z,y) 


z„=rL—y£ + i_ n 
Vn=>-L—ze- |_i_ n 


£-1 £ 

[I — n+1 Z£—n) n d(zi — n+1 Vi— ra+l) 

n=l n=l 

£-1 £ 

= n 9{yi+i - Zi) 0(zi - yi) 

i— 1 2—1 


®e{z,y). (7.6) 


It is also easy to see that 
£-1 2 

n(E ^ s (L-yn + i)^s(L-z n )) 

n =1 s=l 


Zn=>L—y£ + 1 _„ 
yn=>L — Z(i + 1 _„ 


£-1 2 

n(E ^i(^£-n)^s(2/£ —n+l)^ 

n=l s=l 

£-1 2 

=n(E ^(^^(yn+i)). (7.7) 

71—1 S=1 
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Thus, we arrive at 


rL £ 

4>{( T 2i) jk {u)) = / n \e lu<yZn ~ Vn) dz n dy n j @e(z,y) ■ (yi)^j(ze) 

•'° n=l 1 


t-1 2 

*n(£ Tt(z n )^( yri+1 )) := {Tu) kj (u) . (7.8) 

n —1 s=l 


Similarly, using the explicit representations for other operators Tjk(u) one can prove that (17.21) 
implies (17.11) . 


8 Zero modes 


A method of calculating form factors of local operators in GL(3)-invariant models was developed 
in 0- This method is based on the use of partial zero modes of the monodromy matrix entries 
Tij{u) in the composite model consisting of two partial monodromy matrices (12.121) . In spite of 
this approach can be applied to a wide class of integrable models, it should be slightly modified 
in the case of the TCBG model. The matter is that it was assumed in [12J] that the monodromy 
matrix T{u) goes to the identity operator at |u| —>• oo. This restriction is not very important, 
however it leads to minor changes in the case of the TCBG model. 

Observe that a monodromy matrix T^ a \u) constructed by the L-operator m possesses 
the property mentioned above. Indeed, we can define local T-operators Ln\u), (n = 1,... ,N) 
by equations (HUD and QT2J) , where the operators (i^ and are respectively replaced with t//,. (n ) 
and ci\,(n) with the commutation relations [ai(n), a\.(m)] = 5 nm 5ik- Then we can set 


T^ a \u) = u~ n L^\u) 


4 a) («), 


( 8 . 1 ) 


and this matrix obviously has an asymptotic expansion 

T {a \u) =I+-T (a) [0] + O(u" 2 ), u —^ oo. (8.2) 

u 

Therefore we can define zero modes of this monodromy matrix in a standard way 

T (a) [0]= lim -(T (a) (u)-I). (8.3) 

u—> oo C v 7 

However, passing from L-operator (14.11) to L-operator (14.91) we have multiplied L^ (it) by 
the matrix J = diag(l, 1,-1) (see (14.51) 1. This led to the fact that the monodromy matrix T(u) 
(14.81) in the continuous limit has essential singularity at infinity. Therefore, in the case of the 
TCBG model, the definition of zero modes needs to be clarified. We do it in this section and 
consider an asymptotic expansion of the monodromy matrix entries (it) at large value of the 
argument. For this purpose we use the integral representations for Tjj(it) obtained in section [6j 
If it —>• oo, then the expansion for the monodromy matrix contains multiple integrals of 
quickly oscillating exponents. Methods of calculating quickly oscillating integrals are well known 
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(see e.g. [36|, [37[). In our case the integration domain of every integration variable is a finite 
interval [0, L], therefore one of the simplest ways to obtain the asymptotic expansion of Tij(u ) is 
the integration by parts. Using this method one can easily show that single and double integrals 
give 1/rt-behavior, while all the terms with l > 1 give contributions of order o(u~ 1 ). Therefore, 
in order to find zero modes it is enough to take only the first nontrivial terms of the expansion 
for T(u). Then we have 

,-L 

Tij(u) = Sij + x / e iu ^9(z - y)*\{z)*j(y) dzdy + 0(x 2 ), i,j = 1,2, (8.4) 

Jo 

T 33 (u ) = e iuL + xe iuL ( L e iu ^9{z - y) (^(y)*^) + tft(y)tf 2 (z)) dzdy + 0(x 2 ). (8.5) 

Jo 

T a (u) = -iyfr /V^J(y)dy + 0(x 3 / 2 ), i = 1,2, 

Jo 

i(u) = i^xe mL [ e~ my ^j(y) dy + 0(x 3/2 ), j = 1 , 2 , 

Jo 


T 


3 j\ 


( 8 . 6 ) 


(8.7) 


All the terms denoted by 0(x 2 ) or 0(x 3 / 2 ) give contributions 0(u 2 ) as u —>• oo, and therefore 
they are not important. Integrating by parts we obtain 


L y 


(u) = + — f ^\(y)^j(y) dy + 0(i 

u Jo 


— 2 \ 


i,j = 1,2, 


T 33 (u) = e iuL - — e" 
u 


[ {V\(y)*i(y) + *l(y)* 2 (y)) dy + 0{u~ 2 ). 
Jo 


Ta(u) = _^( e »Upt (L ) _ vj/T (0) ) + 0(u- 2 ), 


T 3j (u) = (^(L) - e iuL tyj(0)) + 0(u ~ 2 ), 


Now we define zero modes as follows: 


Tij[0] = lim -(Tij(u) - 8 ^) = - / (y)^j(y)dy, 

u^-oo C J o 


* = 1 , 2 , 

3 = 1 , 2 , 

hJ = 1 , 2 , 


( 8 . 8 ) 

(8.9) 

( 8 . 10 ) 

( 8 . 11 ) 

( 8 . 12 ) 


(recall that c = —ix). This is the same definition as for the models considered in [12}]. The zero 
mode T 33 [ 0 ] is slightly differently: 


X33[ 0 ] = lim - (e luL T 33 {u) - l) = f (^l(y)^i(y) + Tj(y)T 2 (y)) dy, 

u^oo C J 0 

and thus, T n [0] + T 22 [0] = -T 33 [0]. 


(8.13) 
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Looking at (I8.1()D . (18.111) we see that actually we have two types of zero modes for these 
operators. We call them left and right zero modes and denote respectively by T^'- } [0] and [0]. 
Then 




As = 1,2, 


(8.14) 


and 


l\ X 


= JSco7 r « ( “) = iJZ *' (£) ’ 


3 = 1=2. 


(8.15) 


The sums [0] + [0] play the same role as the zero modes of the monodromy matrix of 

the type (18.11) . (|8.2I) . It is known, in particular lij, 3f|, that some of zero modes Tij[ 0] annihilate 
on-shell Bethe vectors: 


i( a ) 


i,b{u,v) = 0, i>j. 


Similarly one can check that 


(4 L) [0] +T 3i 


(r) 


,,b( u ’V) = 0, j + 3, 


(8.16) 


(8.17) 


provided B a b(u, v ) is an on-shell vector. In order to prove (|8.17|) it is sufficient to use the formulas 
of the action of Tij{u ) onto Bethe vectors [25] and to consider there the limits u —>• ±ioo like in 

(EED. 

Finally, the obtained formulas for the zero modes allow us to study form factors of the local 
operators in the framework of the composite model (12.121) . Indeed, let in ()2.12|) the partial 
monodromy matrix T^\u) corresponds to an interval [0,x], where x is a fixed point of the 
interval [0, L\. Then the partial zero modes T^[0] and T-- : [0] are given by (18.121) (18.151) . 
where one should replace everywhere L by x. In particular, we obtain 


*!(*)*:; (®) = = 7Z-JZ }™L u ( T ij\ u ) - <%)> ij = h 2, 


dx 


ix dx u->o o 


$»=ViT^[0] = 


( 1;R )rm — 1 n, 


lim u T 3 L ; ( u ), 


X u—}+ioo 


1 


®fc(z) = ^ T k 3 iH; [°] = lim e T kl \ U )> 


K u^r—ioo 


3 = 1 , 2 , 

k = 1,2. 


(8.18) 


Thus, the problem of calculating the form factors of the local operators in the TCBG model is 
reduced to the evaluating the form factors of the partial zero modes ^ [0] and T^ 1,R) [0]. 


Conclusion 

In this paper we gave a description of the TCBG model in the framework of the algebraic Bethe 
ansatz. The main goal was to prove that the lattice L-operator (14.91) correctly describes the 
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TCBG model in the continuous limit and allows one to find the zero modes of the monodromy 
matrix T(u). This goal is successfully achieved. In order to calculate the form factors of the 
fields 'l'j(x), T|(x), and their combinations \L,}(x)\Lj(x) we can use now the method of 12]. 
Actually, a part of results can be predicted already now. Indeed, the definition (18.121) of the 
zero modes T,j [0] for i, j = 1,2 coincides with the definition used in [12]. Therefore the form 
factors of the operators \E'j(a;)\E f . ; (x) in fact are already computed. The calculation of the form 
factors of the fields 4 l(x) and 'l't(x) should be slightly modified. However, in this case the 
modification affects only the limit u —> oo, but does not affect determinant representations for 
the partial zero modes. We consider this question in details in our further publication. 
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